Current driven magnetization dynamics in helical spin density waves 
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A mechanism is proposed for manipulating the magnetic state of a heUcal spin density wave using 
a current. In this paper, we show that a current through a bulk system with a helical spin density 
wave induces a spin transfer torque, giving rise to a rotation of the order parameter. The use of spin 
transfer torque to manipulate the magnetization in bulk systems does not suffer from the obstacles 
seen for magnetization reversal using interface spin transfer torque in multilayered systems. We 
demonstrate the effect by a quantitative calculation of the current induced magnetization dynamics 
of Erbium. Finally we propose a setup for experimental verification. 
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The possibility to manipulate the magnetization of a 
magnetic material using the spin transfer torque (STT) 
was proposed from theoretical considerations ten years 
ago by Slonczewski El and Berger 0. Since then STT 
has been a very active area of research boosted by its 
potential technological application in magnetic random 
access memories. |3| The phenomenon has been studied 
in layered materials where it is considered to be an in- 
terface effect occurring when a current passes an inter- 
face between two regions with different magnetization 
direction.^] Experimentally, it has been shown that the 
current induced torque is able to switch the magnetiza- 
tion of a ferromagnetic layer in a nano pillar containing 
a ferromagnctic/non magnetic/ferromagnetic trilayer.^Sj 
An obstacle for using the STT in technological applica- 
tions is the magnitude of the smallest current needed to 
switch the magnetization. Theoretical models of trilayer 
systems have shown that the size of this critical current 
is determined by the magnetic anisotropies.j^ 

In this work we consider materials with a helical spin 
density wave (SDW), a magnetic ordering which math- 
ematically can be described as a spin spiral (SS).0, 
In a SS the direction of the magnetization rotates as one 
moves along the SS wave vector. Without loss of gener- 
ality we consider the magnetization rotation axis to be 
parallel to the SS wave vector and will refer to it as the 
spin spiral axis. In comparison with a multilayer a SS 
can be considered as a self assembled magnetic multi- 
layer system where the atomic layers perpendicular to 
the spiral axis are ferromagnetic, but where the magneti- 
zation of the different layers are non parallel with respect 
to each other. We predict that a new kind of STT will be 
induced by a current parallel with the spiral axis. This 
phenomenon is a bulk effect in contrast to the STT in 
a multilayered nano pillar. An advantage with a bulk 
system is the possibility to better control the magnetic 
anisotropics that interfere with the STT in an artificial 
multilayer. 

First we give a general argument for the occurrence of 
STT in SS systems and then we estabhsh the effect by 
a quantitative calculation of the current induced STT in 



Er. Er is chosen since it is one of the many well known 
examples of a helical SDW in the rare earth (RE) se- 
ries (see e.g. Ref. 0). The calculation is performed by 
combining first principles calculations and semiclassical 
linear response theory. 

The STT in multilayered systems can be understood 
by considering a trilayer system A/B/C of two ferromag- 
netic layers A and C, separated by a non magnetic layer 

B. If the two ferromagnetic layers A and C have non par- 
allel magnetization, a current passing perpendicular to 
the layers, from layer A into layer C will change its spin 
polarization when passing into layer C.|4i] The change in 
spin polarization will cause a non zero spin flux into layer 

C, i.e. the spin current entering layer C carries a differ- 
ent spin than the spin current exiting layer C. The net 
spin flux into layer C acts as a torque on the magnetiza- 
tion of layer C. In a SS system a current along the spiral 
axis passes through atomic thick layers with non paral- 
lel magnetization. The spin of the electronic states that 
carry the current rotates coherently with the local mag- 
netization of the spiral, The precessing motion of the 
spins induces a torque on the magnetization that tends 
to rotate the spin spiral. In the ground state the net STT 
from all the electronic states cancels out since there is an 
equal number of electronic states with a current parallel 
and anti parallel to the spiral axis. But if a net current 
flows along the spiral axis no cancellation will occur and 
a STT will act on the total angular momentum of the 
atoms in the material. The torque will cause the mag- 
netization to rotate with a frequency proportional to the 
applied current, or equivalently make the SS slide along 
the spiral axis. 

In order to quantify this effect we calculate the STT 
from a method similar to the method used in Ref. 
for multilayer systems. The calculations are based on 
calculations of the spin flux using the spin current density 
tensor Q (and the STT is investigated both for model 
systems and by using first principles calculations of real 
systems). Q is given by 



Q„k(r) = TrRe [ijl^{r) S ® v^„k(r)} 



(1) 
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where S is the spin operator, v is the velocity operator 
and '0„k are the spinor wave functions with band and 
wave vector indexes n and k respectively. The torque 
exerted by an electronic state with band index n and 
wave vector k on the angular momentum within a volume 
V enclosed by the surface S is given by its spin flux into 
V, 

/^Q„...S^-^. (2) 

In the above equation the torque is expressed as the 
change of total angular momentum 9J„k/9t, within vol- 
ume V. 

The volume of interest is arbitrary, for the STT in 
multilayer systems one considers the torque acting on a 
magnetic layer. For SS systems the STT on individual 
atoms can be calculated using non overlapping spheres 
centered at the atoms. The total torque on an atom is 
obtained by summing the torque from all the occupied 
electron states. If the system is in its equilibrium ground 
state the sum of all torques will be zero due to time 
invariance. However, if an external electric field E is 
applied the occupation of the states at the Fermi surface 
(FS) will change and the total torque on an atom can be 
calculated from a FS integral derived from semi classical 
Boltzmann linear response theory, 

dJ _ Vc re^ f dJ^ dSnk , . 

dt [2ny h ^ Jfs cit |Vke„k| 

where Vc, t, e and e„k are the volume of the system, 
electron relaxation time, electron charge, and band ener- 
gies respectively. The sum is over all bands crossing the 
Fermi level. The above equation defines a linear relation 
between the torque and the external field, 

f^xEA"E. (4) 

n 

A similar expression is obtained for the resistivity, relat- 
ing the current density with the external field. 

j =r^B„E (5) 

n 

— 1 re -k ^ /" v7 -ciN dSnk 

i^^) ^Jfs |Vke„k| 

Combining these two equations gives a linear relation be- 
tween the torque and the current density where the un- 
known electron relaxation time t has been canceled, 

a^ = (EAn)(EBm)"'j = Cj. (6) 

n m 



The above equation for the torque current tensor C will 
now be evaluated for a real SS system, the helical SDW 
in the RE metal Er. Erbium has a complex non-coUinear 
magnetic structure which is strongly temperature depen- 
dent. Bulk Er has an hep structure with c=5.585 A and 
a,=6. 56 A. [13] There is a rich variety of non-collinear mag- 
netic structures and ordering vectors over different tem- 
perature ranges. Below 20K it has a conical SS, between 
53.5 and 85K a longitudinal SDW and between 20K and 
53. 5K is there an intermediate magnetic structure (see 
Ref. m. 

Although the formalism is valid for conical SS, in this 
work we will focus on planar spin spirals. All the mate- 
rial specific quantities of Er used in the calculation of the 
matrix C were calculated from first principle density func- 
tional theory. The calculation of Er was made using the 
full-potential augmented plane wave plus local orbitals 
(FP-APW+lo) method as described in Ref. [H The lo- 
cal spin density approximation (LSDA) as parametrized 
by von Barth and Hcdin was used without use of any 
shape-approximation to the non-collinear magnetization, 
i.e. charge and magnetization densities as well as their 
conjugate potentials are allowed to vary freely in space 
both regarding magnitude and direction. A set of 248 k- 
points was used for converging the electron density. The 
SS was treated using the generalised Bloch theorem. The 
4f-electrons were treated as core electrons. A large set 
of SS wave vectors q along he out-of-plane axis in the 
hep lattice were calculated and an energy minimum was 
found for q~0.20 2tt/c. The results are in agreement with 
previous calculations made for helical SDW in RE.0 

The spin current density tensor was calculated at the 
surface of the augmentation spheres of the atoms where 
the APW expansion can be written as a sum of plane 
waves. 



G 

where a and (3 are the up and down spinors respectively 
and G are the reciprocal lattice vectors. The plane wave 
coefficients a and b are obtained from the first principles 
calculation. The spin flux into a sphere with radius R 
centered at an atom at cite r„ is for plane waves given 
by the expression 
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Qnk • dS 
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a^sa a;k,Ga«k,G'e-'(°-°')"^"ji(|G - G'\R) (G' + k - q/2).(G^G0 
ats/3 a:k,Gfenk,G'e-^(°-°'-'i'"-"ji(|G - G' - q\R) (G' + k + q/2).(G -'G^ - q) 
/3tsa 6:k,G«nk.G'e-^(°-°'+'i)'-"ji(|G - G' + q\R) (G' + k - q/2).(G -'G^ + q) 
/?ts/3 6;k,Gfenk,G'e"''''"°')'""Ji(|G - G'\R) (G' + k + q/2).(G^G' 



(8) 



where m is the electron mass and ji is the first spherical 
Bessel function. In order to get a correct description of 
the FS used in Eqs. © and © a 41 x 41 x 41 k-point 
mesh was used to cover the first Brillouin zone. 

The torque current matrix C is evaluated for an Er 
atom situated at a site with magnetization direction 
[100] and SS wave vector [OOq], 




C = ?i I 0.5 




From the structure of the C matrix we conclude that 
the torque induced by a current along the spiral axis 
causes the SS to translate along the spiral axis, which 
is equivalent to a rigid rotation of the spiral. In order 
to calculate the rotation frequency as a function of the 
current density we also need to consider the total angular 
momentum of the Er atoms J — 15/2 given by Hund's 
rules as the sum of orbital and spin angular momentum 
L = 6 and S* = 3/2 . The torque calculated using the 
C matrix will cause the total angular momentum of the 
atom as well as the SS to precess with 0.07 GHz if a 
current of lO^A/cm^ fiows along the spiral axis of bulk 
Er. 

By analyzing the quantities in Eq. JBJ one finds that 
the main contribution to the C matrix comes from the 
the band whose FS is shown in Fig. ^ This surface is 
the remains of the FS which has the nesting features that 
drives some of the RE solids to have a helical SDW. 9] 
The large contribution from this surface is due to that 
it has a large area being perpendicular to the spiral axis 
that contributes to the electron transport along the axis. 
An estimate of the STT induced by the electronic states 
on this FS can be obtained from the following consider- 
ations. As the conduction electrons flow along the spiral 
axis the component of their spins that are perpendicu- 
lar to the SS axis i.e. for planar SS parallel to the local 
magnetization direction will be precessing around the SS 
axis. 10 The size of the parallel component P can be 
estimated by the expression 




P(k) = 



m(r) ■ Sk(r) 
|sk(r)| 



(9) 



FIG. 1: The Fermi surface of Er that contributes the most to 
electron transport along the SS axis. The spiral axis is along 
the z direction. 



with 

Sk(r) = vlWsV'k(r), 
m(r) = [cos(qr), sm(qr), 0], 

where fn(r) is the local magnetization and (...) means 
space average. P is a measure of the local spin polariza- 
tion and the value of P for the states at the FS in Fig. 
n is shown in Fig. [21 In Fig. [21 is the local spin polar- 
ization P calculated using the interstitial region between 
two atomic planes i.e. the space average in Eq. ^ is 
not done over the whole unit cell. From Fig. [21 we es- 
timate that the spin of the electron states at the FS on 
average are tilted 30° from the spiral axis in the oppo- 
site direction of the local magnetization since the average 
spin polarization P for these states is —0.5. The current 
carried by an electronic state transfers spin upon pass- 
ing an atomic layer perpendicular to the SS axis. This 
is since its parallel spin component will rotate with qir 
[rad] for each layer the current passes. For states on the 
FS a spin of {h/2)'KPq [Js] will be transfered per layer. 
A current of l[A/m^] along the spiral axis induces a STT 
of {h/2)2TTPqA/e [J] per unit cell, where e is the electron 
charge and A = a'^V3/2 is the cell area. The STT in- 
duced by lA/m^ along the spiral axis causes the SS to 
rotate with (PqA)/ (iJe) [Hz]. This estimate of the rota- 
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FIG. 2: Projection of the Fermi surface of Fig. Q on the x-y 
plane. The spin polarization parallel with the local magneti- 
zation direction P is given by the colour code, 1 means that 
the spin polarization is parallel with the local magnetization 
direction, for details see text. 



tion frequency gives four times the value obtained from 
the C matrix but catches the order of the effect. 

For STT in trilayer systems there is a critical cur- 
rent density required to switch the magnetization which 
magnitude is governed by magnetic anisotropies. 
For SS systems there will also be a critical current 
needed to overcome the energy barrier due to magnetic 
anisotropies. This implies that the spiral rotation fre- 
quency will scale linearly with the current along the spi- 
ral axis for currents larger than some critical current of 
the system. But since the STT in SS systems is a bulk 
phenomenon the critical current could be smaller then 
the critical current for trilayer systems. Moreover, for 
the RE systems the anisotropies can be manipulated to 
a large degree, by varying the RE species by alloying, 
while keeping the helical SDW magnetic order. \vi\ 

We suggest that the rotation frequency of the SS as a 
function of current density can be experimentally veri- 
fied in the following way. Consider a nano pillar with a 
layer of a nonmagnetic material in between a SS layer, 
with its spiral axis perpendicular to the interface, and 
a ferromagnetic material as shown in Fig. |31 In princi- 
ple all these three layers can be RE based, e.g. with Gd 
as the ferromagnet and Lu as the non-magnetic species. 
This system will in similarity with GMR and TMR sys- 
tems have a conductance perpendicular to the layers de- 
pendent on the relative magnetization directions of the 
interfaces. 01 If the magnetization of the ferromag- 




FIG. 3: Device for experimental verification of STT in a 
material with a helical SDW. 



netic material is fixed and the magnetization of the SS 
layer precesses then the conductivity perpendicular to the 
layers will be time dependent. The oscillating conductiv- 
ity can be used as a probe of the rotation frequency of 
the SS. 
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